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Original expressions for heat capacity CV and its components, vibrational and configurational components of 


thermal expansion coefficient 


helium 4He were calculated. 


were established. The values of CV, Cvib, Cconf, avib and aconf for water and 
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Introduction 


The heat capacity at constant volume (Cy = dU/dT) 
according to the definition of internal energy (U = Uxint+ 
Upot) has two components. In expression [1] Cy and ap 
components associated with various types of motion of 
the molecules and their component parts (vibrations, li- 
brations, translations, etc.) are called «vibrational» (Cvp, 
Qyib) , and the components associated with molecular 
configuration changes — «configurational» (Cont, Ocont): 

Cy=Cyin +Coone (1) 
Ap=Avib + Ocont (2) 

Phenomenological thermodynamics relates heat ca- 
pacity at constant pressure Cp and constant volume Cy by 
ratio 


Cp - Cy = — (3) 
YT 


where VY) — molar volume, T — temperature (K), a) = 
(1/V)-(6V/0T)p — isobaric thermal expansion coefficient, 
yr = -(1/V)-(ƏV/ðP)r — isothermal compressibility factor, 
(n= TARV / Yr ) — heat capacity caused by thermal 


expansion. 
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Activation of mechanical degrees of freedom taking 
place in the substance under the given conditions (in- 
cluding molecule translations) increases the volume, so 
«vibrational» thermal expansion coefficient should be 
positive (Qip> 0). 

Expression (3) with (2): 

Aino rs T Aono + 2T Ay» AcontM0 
YT YT YT 

In order to determine the type of heat capacity com- 
ponents we consider the following: 

- pand Qon¢ Should have opposite signs, as there 
are states of matter (water, helium, ...) in which œ= 0; 

- expression for Cp and Cy should include aj, and 


Cp oo 


(4) 


Qconfs 

- component (Ta, / yr ) should go to the expres- 
sion for Cy (increase in oscillation amplitude leads to 
increase in the kinetic component of internal energy) and 
C; 
volume is observed); 

- expression for Cy should include negative sign, as 
value Cy> 0, and Qyib, Gone have different signs. 

Given these conditions and expressions (1) - (4) we 
get: 


nem = Tapn [yr (thermal expansion of the system 
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_ TAM m T aconto 


Cp (5) 
YT TE 
Cy= = 2TVoQyib cont VT (6) 
Ta? V, 
Coonf T “cont : (7) 
YT 
TaVo TaVo 
C= — t (8) 
YT YT 


Replacing sont with (œp — Qi») in (5) we obtain a 
quadratic equation the physical meaning of which has the 
following solution: 


Ain = (9) 


Table 1 Vo, Cp,yr and a of water at atmospheric pressure [3-5] 
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Expressions (1) — (9) allow us to find value Cy, Cyin, 
Coons Aib ANd Qeons, if we know values Vo, Cp, yr and ap. 

Calculated values Qyij, Qconts Ceonf and Cy, for water 
H2O and liquid helium “He are shown in Tables (1-5) and 
Figs. (1-4). For this research, condensed state of water 
and helium are interesting because dependence Vo= fT) 
and dependence œ= f(T) of these substances have ex- 
treme points. 


Water 


Reference values Cp, Cy,Vo, Œp , yr [3-5] and calcu- 
lated values Qcont, Cy, Cvibs Coont and Crnerm for protium 
water from 273.15 to 373.15 K are shown in Tables 1 and 2. 

From Tables 1 and 2 we can see that the vibrational 
and configurational components of thermal expansion 


a, 107K" 
T,K Vo, 10° m°/ mol Cp, J/ mol-K 1,107! Pa! 
Op Qvib Qeonf 
273.15 18.0184 76.028 50.886 -0.681 19.482 -20.163 
277.13 18.0165 75.850 49.514 0.000 19.393 -19.393 
283.15 18.0214 75.595 47.811 0.879 19.253 -18.374 
293.15 18.0484 75.397 45.895 2.066 19.087 -17.021 
303.15 18.0948 75.307 44.774 3.031 18.981 -15.950 
313.15 18.1574 75.289 44.243 3.890 18.949 -15.059 
323.15 18.2341 75.307 44.260 4.514 18.922 -14.408 
333.15 18.3236 75.361 44.500 5.290 19.006 -13.716 
343.15 18.4254 75.451 45.210 5.830 19.078 -13.248 
353.15 18.5388 75,595 46.160 6.472 19.237 -12.765 
363.15 18.6634 75.137. 47.310 6.990 19.375 -12.385 
372.85 18.8010 75.950 50.404 T332 19.857 -12.325 
373.15 18.7990 75.973 48.900 T3537 19.599 -12.062 
Table 2 Heat capacity of water and its components at atmospheric pressure 
C, J/ mol-K 
T,K 
Cy, ref Cy, cale Cyb Cront (Tho a yr 
273.15 76.028 75.983 36.664 39.320 0.045 
277.13 - 75.850 37.925 37.925 0.000 
283.15 75.433 75.513 39.480 36.033 0.082 
293.15 74.820 74.905 41.506 33.399 0.492 
303.15 74.226 74.181 43.014 31.168 1.126 
313.15 73.469 73.344 44.201 29.144 1.945 
323.15 72.640 72.594 44.956 27.638 2.713 
333.15 71.740 71.522 45.715 25.808 3.839 
343.15 70.320 70.698 46.151 24.547 4.753 
353.15 69.848 69.654 46.544 23.110 5.941 
363.15 68.875 68.757 46.781 21.976 7.000 
372.78 - 68.062 46.939 21.123 7.888 
373.15 67.884 67.825 46.954 20.871 8.149 
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have opposite signs; dependence a,;,= AT) has a mini- 
mum value in a temperature range of (30-50)°C ap- 
proximately where the value of heat capacity Cp is mi- 
nimal It follows from Table 1 that the increase in ther- 
mal expansion coefficient of water in the range of (Tine — 
T») is caused by a decrease in modulus Qone. 

According to Morgan and Warren [6], melting ice and 
heating water from 0 to 62°C simultaneously increase the 
number of nearest neighbors of the molecule and the av- 
erage intermolecular distance. Decrease in the volume 
of condensed water during melting and heating from 0°C 
to 3.98°C is explained by the fact module Qeon¢ in these 
processes is greater than aj, [2] (volume reduction due 
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to increase in the number of nearest neighbors molecule 
is greater than increase in volume because of the growth 
of intermolecular distance). At a temperature of maxi- 
mum density (3.98°C) water |Qeond=Qyip. With further 
heating, the volume of water increases as Qyip > |Qcont- 
Reference value a, and calculated values Oyi , cont 
in liquid water at saturation line from 373.15 K to sub- 
critical T>0.96 T, (T.— critical temperature) are shown in 
Table 3. Type of dependences a= f(T) and @con AT) 
are shown in Fig. 1. The values ap and yr were deter- 


mined using (Z) 5 (Z) [3-5]. 
OT Jp OP Jr 


Table 3 Reference mp and calculated values Cyib, Coons, Ovib and Ocone at saturation line in liquid water 


T, (K) Civ, J/ mol-K Ceont, J/ mol-K ap, 107 K” Qyib, 104 K! Acon 104 K! 
373.15 46.954 20.871 7.537 19.599 -12.062 
393.15 47.257 20.403 7.900 19.904 -12.004 
403.15 47.314 19.085 8.650 20.346 -11.696 
413.15 47.286 17.964 9.350 20.806 -11.456 
423.15 47.397 17.529 9.750 21.194 -11.444 
433.15 47.391 16.751 10.400 21.772 -11.372 
443.15 47.595 16.540 10.800 22.304 -11.504 
453.03 47.047 14.810 12.151 23.315 -11.164 
485.52 46.755 12.740 14.876 26.019 -11.143 
506.99 46.864 11.587 17.266 28.609 -11.343 
523.48 46.704 10.613 19.645 31.199 -11.554 
537.07 46.683 9.888 22.089 33.961 -11.872 
548.71 46.796 9.279 24.720 36.947 -12.227 
558.95 46.795 8.683 27.603 40.180 -12.577 
567.95 46.734 8.086 30.848 43.756 -12.908 
576.46 46.908 7.604 34.492 47.839 -13.347 
584.11 47.356 7.194 38.689 52.550 -13.861 
591.19 47.663 6.769 43.529 57.938 -14.409 
597.79 48.084 6.358 49.320 64.351 -15.031 
603.96 48.779 5.989 56.355 72.130 -15.775 
609.78 49.619 5.634 65.200 81.903 -16.703 
615.27 51.427 5.490 77.200 95.654 -18.454 
617.90 51.112 5.079 84.900 103.634 -18.734 
620.47 53.193 5.177 93.000 113.055 -20.055 
622.97 53.086 4.723 104.100 124.435 -20.335 
625.41 53.758 4.507 115.900 137.114 -21.214 
627.79 60.559 5.284 129.900 154.735 -24.835 
630.11 58.086 4.309 152.300 176.704 -24.404 
632.37 57.464 3.708 180.500 205.401 -24.901 
634.59 58.027 3.251 220.000 246.115 -26.115 
638.86 55.811 1.963 377.000 404.491 -27.491 


128 
80 
60 
Oni 
40 
20 
0 T T T T T T T IT. K 
240 300 350 40) 50 500 550 600 650’ 
ao | ane et 
Qeont 
—40 
Op, LO“K" 


Fig. 1 Temperature dependence of Calculated values a,j, and 
Ocone Of water at saturation line 


Helium 


For calculation of values Qyib, Gon Coon and Cyip for 
liquid helium “He we use Cp, Cy Vo, œ and yr (Tables 4-6) 
[7-9]. Thermal expansion coefficient according to [7, 8] 
in the temperature range of (1.15-1.20) K (Fig. 2) and in 
the temperature range of (2.15-2.20) K is equal to zero. 
It allows us to define vibrational and configurational in- 
puts to the thermal expansion coefficient of liquid “He. 

Liquid “He can exist in two fractions: He II and He I. 
Helium II — a fraction with a superfluid component (su- 
perfluid), the proportion of which varies from 1 for 0 K 
to 0 for 7,=2.1768 K (the temperature of 4- transition) 
[8]. 
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Fig. 2 Temperature dependence of thermal expansion coef- 
ficient 


Tables 4 and 5 provide reference and calculated values 
of liquid “He II at atmospheric pressure. 

The values of heat capacities Cp and Cy, Cy and ap 
components of liquid helium at saturation line are given 
in Table 6. Fig. 3 shows the dependence Otcon¢ =f (T) and 
Fig. 4 - dependence of configurational and vibrational 
heat capacity of liquid helium at saturation line. 


Results and discussion 


The separation of heat capacity Cy and thermal expan- 
sion coefficient ap into components is good for studying 
the mechanism of formation of thermal and structural 
properties of the liquid state. For example, Tables 1-3 
show that value Qon¢ decreases for water being heated 
from melting temperature (Tnmen) to Tmax¥0.7857,508 K. 


Table 4 Vo, ap,y7, Cp, and Cy of liquid tHe II at atmospheric pressure [8] 


T,K Vo, 10° m°/mol ap, 104 K" Yr, 10° Pa” Cp, J/ mol-K Cy, J/ mol-K 
0.6 27.580 1.749 12.147 0.017733 0.017729 
0.8 27.582 3.213 12.156 0.084307 0.084286 
1.0 27.583 3.294 12.193 0.42047 0.42043 
1.05 27.584 2.495 12.210 0.58638 0.58638 
1.10 27.584 1.225 12.225 0.79320 0.79316 
1.15 27.584 0.04805 12.249 1.04364 1.04364 
1.20 27.584 -2.215 12.264 1.34183 1.34179 
1.25 27.584 -4.523 12.290 1.68626 1.68622 
1.30 27.583 -7.342 12.313 2.10721 2.10705 
1.35 27.582 -11.17 12.335 2.63215 2.63179 
1.4 27.580 -13.87 12.374 3.23962 3.23902 
1.6 27.566 -35.79 12.492 6.21083 6.20643 
1.8 27.539 -65.82 12.641 11.73602 11.71881 
1.9 27.516 -87.32 12.720 15.48686 15.45564 
2.0 27.488 -112.8 12.943 20.82753 20.77349 
2.05 27.471 -141.6 13.169 23.70460 23.61894 
2.10 27.449 -177.0 13.267 26.72456 26.58847 
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Table 5 Vibrational and configurational inputs to Cy and ap of liquid “He I 
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T,K vib LOW K! Acong 107 K” Cvi», J/ mol-K Coons, J/ mol-K 

0.6 81.545 -79.796 0.009054 0.008674 

0.8 153.987 -150.774 0.043023 0.041265 

1.0 154.014 -150.720 0.053651 0.051389 

1.15 448.907 -448.859 0.52188 0.52176 

1.20 497.466 -499.685 0.66788 0.67390 

1.40 713.519 -727.389 1.58803 1.65099 

1.60 919.775 -955.565 2.98241 3.22391 

1.80 1189.93 -1255.75 5.53515 6.18366 

2.00 1508.38 -1621.18 9.6100 11.1635 

2.10 1662.96 -1839.96 11.87926 14.70921 

Table 6 Characteristics of liquid tHe at saturation line [9] 
TK Heat capacity, J/ mol-K ap, 10° K"! 
Cp, ref Cy, ref Cw C oont Cy, calc Ovib conf 

0.800 0.089738 0.089698 0.04604 0.04366 0.08970 159.593 -155.348 
0.900 0.19877 0.19873 0.1010 0.0977 0.19874 222.524 -218.746 
1.000 0.40106 0.40106 0.2018 0.1993 0.40105 298.081 -296.233 
1.050 0.55706 0.55706 0.2756 0.2751 0.55076 340.084 -339.704 
1.100 0.74128 0.74128 0.3693 0.3720 0.74128 384.757 -386.153 
1.150 0.97944 0.97944 0.4858 0.4936 0.97941 431.860 -435.321 
1.200 1.27123 1.27123 0.6279 0.6432 1.27114 481.131 -486.922 
1.400 3.12083 3.12003 1.52108 1.59881 3.11989 697.426 -714.805 
1.600 6.41617 6.41216 3.09664 3.31590 6.41254 939.191 -971.303 
1.800 11.75964 11.74763 5.61745 6.13142 11.74888 1207.600 -1260.428 
2.000 20.761 20.721 9.690 11.030 20.720 1541.129 -1640.729 
2.100 28.995 28.883 13.1126 15.7704 28.8830 1786.887 -1951.315 
2.150 37.100 36.844 16.1210 20.7235 36.8445 1990.535 -2239.186 
2.200 16.899 16.855 9.015 7.840 16.855 1472.226 -1369.590 
2.400 9.506 9.014 5.112 3.240 9.012 1181.181 -849.306 
2.600 9.286 8.337 5.739 2.599 8.338 1214.368 -757.060 
2.800 9.618 8.125 5.889 2.233 8.122 1286.177 -707.248 
3.000 10.395 8.249 6.215 2.034 8.249 1382.066 -681.733 
3.200 11.592 8.686 6.728 1.958 8.686 1505.551 - 678.676 
3.400 12.940 9.134 7.248 1.888 9.136 1652.886 -683.180 
3.600 14.413 9.510 7.717 1.791 9.508 1826.996 -688.107 
3.800 16.078 9.802 8.133 1.665 9.798 2043.173 -694.489 
4.000 18.088 10.022 8.507 1.515 10.022 2325.440 -703.190 
4.200 20.693 10.199 8.855 1.334 10.199 2716.396 -715.920 
4.400 24.404 10.343 9.197 1.152 10.349 3300.384 -734.475 
4.600 30.380 10.479 9.548 0.933 10.481 4276.202 -760.995 
4.800 42.067 10.627 9.966 0.685 10.651 6247.818 -804.068 


In [10] it is shown that at Tmax volume multiplied by 
pressure of saturated gas (PVo)s, and molecular packing 
coefficient (Ks) in liquid saturation line branch have the 
maxima. Further heating of water (Table 3) increases 
value Qon. Value Adon AT in the temperature range of 
(273.15-373.15) K and at T>0.96T,is considerably high- 


er than in other temperature ranges. This fact shows that 
at low (273.15-373.15) K and high (T> 0.96T,) tempera- 
tures structural changes occur more intensely. In the de- 
pendence Qconr = AT) two distinct temperatures (Tinax® 
0.7857, and Tx0.96T.) may be attributed to temperature 
of quasi-phase transition in liquids [11]. 
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Fig.3 gone =f (T) of liquid ‘He at saturation line 
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Fig. 4 C,» = AT) and Crone = AT) of liquid *He at saturation 


line 


Analysis of dependence Cyonr= AT) (Tables 2 and 3) 
shows that Cone is constantly decreasing with the in- 
crease in temperature of water. Ratio Coon/Cy for wa- 
ter according to our data is 0.518 (273.15 K) and 0.308 
(373.15 K) ), and according to [1] it is approximately 
equal to 0.50 (273.15 K) and 0.35 (373.15 K). 

In [12,13] calculated are the values of radial distri- 
bution of water. It was found that by heating water un- 
der pressure 1000 bar from 298 K to 450-500 K the 
number of nearest neighbors of Z in water molecules 
increases from 4.9 to 5.7, and after further heating - it 
decreases to about 4.3 at 773 K. From Table 3 and Fig. 1 
we see that the value of configurational components of 
thermal expansion coefficient is constantly decreasing 
from melting temperature to 506.99 K. 506.99 K is al- 
most equal to Tmax and is close to the upper limit of 
temperature range of (450-500) K. 

Change of Z may be related to the change of Qeon¢. If 
the increase in Z is accompanied by the decrease of Qons, 
the increase in modulus sonf must be accompanied by 
the decrease of Z. At T>T max modulus Qonf in liquid wa- 
ter is constantly increasing, so we should expect de- 
crease of Z for water molecules in the interval Tmax-Te- 
This assumption may be confirmed by temperature de- 
pendence Ks for water at saturation line [10]: Ks increas- 
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es up to Tmax and then decreases. 

From Tables 5 and 6 we see that modulus Gone for 
liquid *He increases from almost 0 K to T= 4.8 K (Table 
6), except for the interval 7,+3 K. We could expect that 
in He II and He I in the temperature range 3-4.8 K the 
number of nearest neighbors of atoms will decrease (ac- 
cording to [14] the value of Z for “He atoms in the tem- 
perature range of 3-5.1 K decreases). 

Low accuracy in determining the experimental cha- 
racteristics of Helium (Cp, Cy, P, Œp, yr) does not allow 
us to unambiguously interpret the dependence @.on¢ =f (T) 
at high temperatures. For this reason, we used tempera- 
ture 4.8 K as the maximum temperature of liquid “He. 

In the temperature range (0 K - 73) vibrational Cy and 
configurational Coons components of heat capacity in- 
crease; in the temperature range (2.6—4.8) K Cyi, increas- 
es and Cyone for the He I, as it happens to water at 
temperature from 273.15 K to 638.86 K, constantly de- 
creases (Fig. 4). 

In conclusion, the original expressions for heat capac- 
ity Cy and its components were proposed, and the separa- 
tion of thermal characteristics of the substance to «vibra- 
tional» and «configurational» components is good for a 
better understanding of the structural changes in the liq- 
uid state was shown. 
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